Tunable terahertz oscillations in superlattices subject to in-plane magnetic field 
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We present a concept of terahertz oscillations in superlattices generated under conditions appar- 
ently different from standard Bloch oscillations. The oscillations are induced by crossed magnetic 
and electric fields both applied to the superlattice in the in-plane direction. The frequency of these 
oscillations is tunable by the applied fields. 
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I. INTRODUCTION 

The possibility of Bloch oscillations (BOs), i.e. the pe- 
riodic motion of an electron in a periodic system induced 
by a uniform electric field, was mentioned by Zeneai rel- 
atively soon after the basic quantum mechanical theory 
of the solid state was established. 2 Taking a simple pic- 
ture of electrons, which are not subject to tunneling to 
other bands or scattering, Zener predicted an oscillatory 
motion in the real as well as in the reciprocal space with 
the frequency lobo — \e\FA/H, where A and F are the 
system spatial period and the electric field, respectively. 
However, it took a long time before any experimental 
evidence of BOs has been foundi&i The key feature for 
their observation was a pioneering concept of a super- 
lattice (SL) and minibands suggested by Esaki and Tsu£ 
that allowed to overcome problems of a strong electron 
scattering, which made the observation of BOs in bulk 
semiconductors hardly feasible. 

In this paper, we present a simple idea of electron-in- 
plane-oscillations appearing when crossed in-plane mag- 
netic and electric fields are applied to SL. We start with 
a general discussion of properties of SL subject to the 
in-plane magnetic field the effect of which cannot be de- 
scribed within the quasi-classical approximation. The 
findings are illustrated by simple numerical calculations 
based on the standard tight-binding (TB) model&i 

Subsequently, we study influence of an additionally ap- 
plied electric field and conclude that an oscillatory mo- 
tion of electrons is induced. We present two models of 
these oscillations, i) quasi-classical and ii) pure quantum- 
mechanical. The predicted oscillations are compared to 
common BOs. 



II. SUPERLATTICE SUBJECT TO IN-PLANE 
MAGNETIC FIELD 

Let us consider an infinite SL having its growth axis 
oriented along the z-direction, which is described by the 
periodic potential V{z) — V(z + A), where A is the pe- 
riod of SL. The Hamiltonian of an electron in such a sys- 



tem subject to the in-plane magnetic field B = (0, Bn , 0), 
with the vector potential gauge A = (B\\z, 0, 0), reads: 
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+ V(z). (1) 



To find eigenstates of this Hamiltonian, the following 
ansatz for the wave functions is commonly assumed 



XkJ z )- 
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This way the three-dimensional Hamiltonian Q is re- 
duced to one-dimensional H which depends on the pa- 
rameters k r and k„: 
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V(z). (3) 



Hence, when the in-plane magnetic field is applied to the 
SL system, the motion in the x- and z-directions become 
coupled. This coupling can be interpreted as an effect of 
the Lorentzian force. 

To solve the eigenvalue problem defined by the Hamil- 
tonian we can utilize the invariance of this Hamilto- 
nian under transformation H(z, k x ) — > H(z+A, k x — Ko), 
where Ko = \ e \B\\A/h. The whole eigenenergy spectrum 
then reads: 



E n (k x , ky) 



E n (k x ) + ^-, 



(4) 



where E n (k x ) are Landau subbands (n = 1,2,3...) 
which are -Ko~periodic in the momentum k x . 

A straightforward solution of the Schrodinger equa- 
tion including the Hamiltonian @ leads to the functions 

ip^k ( x iUt z ) that do not respect the translation sym- 
metry of the SL along the growth axis. This is caused 
by the additional effective potential parabolic in z, devel- 
oped at finite B\\ in the Hamiltonian Consequently, 
electron momenta k x have to be taken within the interval 
k x G (— oo, +oo) to obtain full set of eigenfunctions J5J. 
Nevertheless, all eigenstates in the n-th subband hav- 
ing momenta k x + IKq (I G Z) are degenerate due to 
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FIG. 1: Electron dispersion curves E n (k x ) calculated in TB 
approximation for three lowest lying subbands n = 1,2,3. 
The parameters used in the calculation are given in the inset. 



the periodicity of E n (k x ) and therefore, the translation 
periodicity can be restored, when an appropriate linear 
combination of these states is taken: 



*k.,k v ,q( x >V> z ) K H e 
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The new continuous quantum number q 6 (— 7r/A, 7r/A) 
thus replaces the discrete index I, Owing to the restored 
translation periodicity, the eigenstates (0) fulfill the con- 
dition: 

V, * " A) = e^+^^Jx, y, z) (6) 

and the considered interval of momenta k x can now be 
reduced to the Brillouin zone k x £ (—Kq/2, Kq/2). Both 
reduced and extended schemes are fully equivalent. 

Hence, the Hamiltonian J3J describes a system with 
periodic dispersions E n (k x ), whose period Kq is tunable 
by the applied magnetic field Bn. This is in contrast 
to a fixed period 2tt/A of the electron dispersion in the 
growth direction of SL at zero Bn. The obtained result 
can be interpreted also as a formation of a 2D lattice 
in the x — z plane induced in SL by a finite magnetic 
field in the y-direction. This lattice has spatial periods 
in x- and z-directions 2it /Kq and A, respectively. The 
magnetic flux trough the unit cell of the lattice is simply 
27r_B|| A/ Kg = h/\e\, i.e. one magnetic flux quantum. 

The form of the Hamiltonian Q implies the appear- 
ance of the periodic dispersion even at a negligible small 
-Bm and thus the standard parabolic dispersion is not at- 
tained in the limit Bn — > 0+. The problem lies in the 
infinite size of the considered SL. When a (realistic) su- 
perlattice with a finite number of wells is taken into ac- 
count, we get only a partially periodic dispersion E n (k x ) 



and the number of minima in this dispersion corresponds 
to the number of wells, see an extreme case of a double 
quantum well in the Appendix. These minima clearly 
disappear at low B\\ and the dispersion approaches the 
expectable parabolic shape. 

Before we further utilize the above drawn conclusions, 
we use a simple TB approximation to get some numerical 
results, which can illustrate the studied problem. Within 
the framework of the TB model, the Hamiltonian @ is 
transformed into the matrix form and readsiiS*ii 

h 2 

Hr.s = ^-(k x + rKofSr.s + t5 r , s±l (r, s £ Z), (7) 
Am 

where the coefficient t (t < 0) characterizes the tunneling 
between adjacent quantum wells (QWs) and where the 
motion in the y-direction is not included, since it is not 
affected by B\\ . Note that this TB Hamiltonian conserves 
the periodicity of the original Hamiltonian @ in momen- 
tum k x . The eigenvalue problem given by the tridiagonal 
Hamiltonian J7J can be easily solved numerically. 

The calculated dispersions for three lowest lying Lan- 
dau subbands E n (k x ) have been plotted into Fig. ^ 
These results demonstrate both the expected periodic- 
ity of subbands in k x and the shape of the dispersion 
curves, which cannot be predicted only from the symme- 
try of the Hamiltonian ©. We see that just one mini- 
mum per interval (k x , k x + Ko) appears. The dotted lines 
in Fig. n show the limit of very weakly coupled wells, 
i.e. t — v 0. The dispersion curves at t = are purely 
parabolic and corresponds to dispersions of electrons in 
isolated QWs. Hence, the widths of individual Landau 
subbands E n (k x ) and the energy gaps between them are 
given by the strength of t and can be also tuned by the 
applied magnetic field Bn . 

Henceforth, we will take account of the lowest lying 
subbands E\{k x ) = E(k x ) only. Such approximation is 
meaningful in strongly coupled SLs, i.e. for high values of 
|i|, when the separation of this lowest subband from the 
higher ones is significant. The energy width Wq of this 
subband can be simply estimated at high magnetic fields. 
We just take W » h 2 {K /2) 2 /(2m) = e 2 B 2 A 2 /8m. Ob- 
viously, this rough approximation fails if Wq « |t|. 



III. SEMI-CLASSICAL MODEL OF 
OSCILLATIONS 

Having the periodic band structure E(k x ) at a given 
fixed magnetic field Bu , we use a semi-classical consid- 
eration to describe the electron motion if an additional 
constant electric field F x is applied in x-direction. As 
the influence of Bn has already been included in the 
discussed energy spectrum, the semiclassical equation of 
motion takes a simple form hk x = —\e\F x and thus k x 
changes linearly in time. Hence, the electron velocity 
v x = hr l dE(k x ) /dk x becomes periodic in time and a 
specific oscillatory motion is generated. This motion is 
schematically shown in Fig. [21 and can be decomposed 
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and thus, taking account of the periodicity in k x , the 
eigenenergies can be easily calculated: 
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FIG. 2: A schematic picture of SL showing an expected tra- 
jectory of the electron under depicted conditions. 



into a steady shift in the SL growth direction and the 
oscillations in the rr-direction. The drift motion in the 
growth-axis direction becomes apparent especially in the 
extended scheme of E(k x ). The used equation of motion 
gives us also possibility to calculate the corresponding 
oscillatory frequency wg,, = 2irF x /(B\\A). Hence, Wg. 
is tunable not only by the electric field, as in the case 
of BOs but by B\\ as well. Both frequencies ujbo and 
Wb\\ are functions of A - but whereas the first frequency 
is linear in A, the latter one has the reciprocal depen- 
dence. The oscillatory frequency wb» can be rewritten 
into u>b« — 2nvd/A, where va = F x /B\\ is the drift veloc- 
ity introduced in 3D for the electron motion perpendicu- 
lar to the crossed electric and magnetic fields^ The ratio 
A jvd is then obviously the time needed by an electron to 
tunnel into the adjacent QW. 

The semiclassical model allows us to determine the 
spatial amplitude of expected oscillations Xq defined in 
Fig. [3 When we make use of the facts that the elec- 
tron position is the time integral of the electron velocity 
and the velocity v x is derivative of the dispersion curve, 
we obtain a simple relation xq = Wo/(2\e\F x ). As the 
subband width Wq varies with .By, the amplitude xq is 
tunable by the magnetic field as well. 



IV. QUANTUM-MECHANICAL MODEL OF 
OSCILLATIONS 



The electron motion in a system with a periodic dis- 
persion E(k x ) = E{k x + Kq) can be treated in a pure 
quantum-mechanical way as reviewed e.g. by Hartmann 
et al. in Ref. The corresponding Hamiltonian is there 
conveniently written in the momentum representation: 



E n = 



1 

KoJ„ 



K 



E{k x )dk x + nfko 



Bu- 



rn 



H(k x ) = E(k x )+i\e\F, 



d 



(8) 



Because the first term of this eigenenergy is constant at 
given Bu , we receive an analog of the common Wannier- 
Stark ladder (n £ Z) discussed in the framework of Bloch 
oscillations i& 



V. REMARKS ON A POSSIBLE REALIZATION 

From a practical point of view, the experiments prov- 
ing emission of the predicted THz radiation can be the 
same as in the case of standard BOs. Both coherent 
and incoherent radiations can be obtained. The coher- 
ent THz radiation induced by BOs is achieved e.g. when 
free electrons arc generated in an undoped SL by a fem- 
tosecond optical pulse ensuring the same phase of all 
electrons The photon laser energy is tuned to an ap- 
propriate electron- hole transition in SL. The same tech- 
nique can be used in our case as well. 

The function of the THz generator could be disrupted, 
if electrons initially localized in the lowest subband tun- 
nel under the effect of F x to the higher subbands. In such 
a case, the one-subband model utilized in both semiclas- 
sical and quantum-mechanical treatments of oscillations 
would not be applicable. We illustrate this obstacle on 
a simple model of a double quantum well in Appendix. 
This model offers the simplest possibility to check the 
intersubband tunneling induced by the electric field. It 
cannot serve as a definite evidence that the tunneling to 
higher subbands is negligible in superlattices, neverthe- 
less, it illustrates that electrons do not noticeably tunnel 
to the higher (antibonding) subband in DQW at Sii un- 
der conditions typical for the THz oscillations predicted 
in SLs. Apparently, further investigations in this direc- 
tion are necessary. 

For a possible realization, we should also check the 
sample design and experimental conditions to observe the 
predicted oscillations in the terahertz region. Assuming 
ujb» ~ ^>bo and the same electric field in both cases, we 
obtain the corresponding magnetic field Bu w h/A 2 \e\ = 
16 T for A = 16 nm. Hence, at Bu < 16 T the oscillatory 
frequency is even higher than for BOs, since cjg^ oc B^ 1 . 
Moreover, having two free parameters F x and Bu we can 
independently optimize cjbu and xq to achieve the maxi- 
mal emitted power. This is impossible for standard BOs, 
since lobo and the corresponding spatial amplitude are 
governed by the applied electric field only. 

The important point in the observation of BOs is the 
achievement of an oscillation period significantly lower 
than is the scattering time due to phonons or plasmons. 
We predict our oscillation for SL systems, where common 
BOs are observed. Therefore, the same or very similar 
damping rates as observed in BO experiments could be 
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FIG. 3: The occupancy \cb\ 2 (\ca\ 2 ) of the bonding (anti- 
bonding) subband in DQW after the reflection of an electron 
on a barrier formed by the lateral in-plane electric field. The 
inset schematically illustrates the coming (IN) and reflected 
(OUT) states in the electron subband structure corresponding 
to the region x < 0. 



expected in our case as well. Hence, the published exper- 
imental evidence of BOsj&iii suggests that the predicted 
B||-controlled oscillations ought to be experimentally ob- 
servable. 

It is interesting to investigate also the direction char- 
acteristics of the expected THz radiation. Since the ra- 
diation is generated by the electron oscillatory motion in 
the x-direction, the radiation should be emitted mainly 
in the plane perpendicular to the cc-axis, i.e. in the plane 
perpendicular to the oscillating dipoles. The predicted 
device can thus be both edge- or surface-emitting. 

An important advantage of presented model in com- 
parison with standard BOs is a fast (in-plane) drain of 
electrons from the structure after they reach the edge of 
SL. This fact can be utilized in a significant enhance- 
ment of the repetition frequency of the generation of the 
coherent THz radiation. 



VI. CONCLUSIONS 

We have investigated behavior of electrons in a super- 
lattice when crossed magnetic and electric fields are ap- 
plied, both in the in-plane direction. We predict a novel 
terahertz oscillations in superlattices that are different 
from Bloch oscillations that appear when the electric field 
is applied in the growth direction of the superlattice. We 
have also found a simple expression for the frequency of 
such oscillations. The suggested realistic design of the 
structures allows preparation of terahertz emitters con- 
trolled by the in-plane magnetic field. 

VII. APPENDIX 

Let us assume a double quantum well (DQW) oriented 
as used in the paper subject to the in-plane magnetic 



field and described within the TB approximation by the 
interwell distance A and the tunnelling coefficient t. The 
thorough theoretical analysis of such a DQW system can 
be found elsewhere&i^ The DQW represents an extreme 
case of a superlattice taken up to now into account. A 
lateral electric field applied in the ^-direction at x > 
forms the potential profile: 

The electron described by the bonding subband wave 
function ty B (k x ,x) with the momentum k x > and en- 
ergy E B (k x ) enters into the system at x < 0. It is re- 
flected by the potential barrier <p(x) and leaves the system 
again at a; < 0. The reflected electrons may be found 
both in bonding states \?b( — k x ,x) and in antibonding 
states ^a{— k xl x) with k' x > being a momentum of the 
corresponding antibonding state. Due to the elastic re- 
flection the energy of antibonding state Ea(— k' x ) equals 
to bonding ones EA{—k' x ) = EB{±.k x ). 

The total wave function ^>(x) of the electron at x < 
thus reads 

, p _ * B (k x ) + c B V B (-k x ) + ca^a(-K) 

yj\ + C B C* B + CAC* A 

where the complex amplitudes of reflected waves c B and 
ca determine the occupancy of respective bonding and 
antibonding subbands. The wave function 'F at x > is 
calculated solving the Schrodinger equation numerically 
and respective c B and ca are established to accomplish 
the damping lim x -> +OQ fy(x) — 0. Results of the calcula- 
tion for different tunnelling coefficients t are shown in Fig. 
13 where model parameters m = 0.067mo, A = 16 nm, 
£>d = 12 T and F x = 1920 V/cm corresponding to the 
oscillator frequency of 1 THz have been used. The en- 
ergy of incoming electrons E B (k x ) well above the mini- 
mum energy of the antibonding subband Ea(0) was used 
to enable the intersubband tunneling. We have ascer- 
tained that the course of \c B \ 2 (\ca\ 2 ) is practically in- 
sensitive to E B (k x ). We find out that increased subband 
splitting strongly damps the intersubband tunnelling and 
\ca\ 2 < 0.01 is obtained at \t\ > 7.5 meV for chosen pa- 
rameters. Analogous results are obtained also for other 
sets of parameters producing oscillations in THz branch. 
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